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Abstract 



In a series of recent publications, a new approach to the non-relativistic reduc- 
tion of the electromagnetic current operator in calculations of electro-nuclear 
reactions has been introduced. In one of these papers Q, the conjecture that 
at energies of a few GeV, the bulk of the relativistic effects comes from the 
current and not from the nuclear dynamics was made, based on the large 
relativistic effects in the transverse-longitudinal response. Here, we explicitly 
compare a fully relativistic, manifestly covariant calculation performed with 
the Gross equation, with a calculation that uses a non-relativistic wave func- 
tion and a fully relativistic current operator. We find very good agreement up 
to missing momenta of 400 MeV/c, thus confirming the previous conjecture. 
We discuss slight deviations in cross sections for higher missing momenta and 
their possible origin, namely p-wave contributions and off-shell effects. 

24.10.Jv, 25.30.Fj, 25.10-Fs 
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I. INTRODUCTION 



Currently, there is a lot of activity in the field of electron scattering from nuclei at 
a few GeV. There exists a broad experimental program of electron scattering studies at 
Jefferson Lab, MAMI, and Bates, aimed at understanding the short range structure of 
nuclei and the properties of nucleons in the nuclear medium. It would be ideal to perform 
theoretical calculations of these processes in a completely microscopic fashion, starting from 
a Lagrangian. Such an approach would contain explicit relativistic treatments of the nuclear 
dynamics and the electromagnetic current, as well as a consistent treatment of the initial 
nuclear state and the in general quite complicated hadronic final state. In practice, for high- 
energy inelastic electron scattering, it is very difficult to calculate both the nuclear ground 
state and the final hadronic scattering state consistently and it is unlikely that a consistent, 
fully microscopic and, therefore, relativistic treatment will be available for medium and 
heavier nuclei in the near future. Currently, there are several approximate calculations 
for particular electro-nuclear reactions available for few-body systems 0,0]. However, it is 
difficult to extend these approaches for the few-body systems to energies above the pion 
emission threshold. With the considerable effort one expects will be put into this field, in 
the next few years calculations for the few GeV regime for reactions with deuteron targets 
will hopefully become available. For heavy nuclei, there are only relativistic mean-field 
calculations available, which do not allow for the investigation of ground-state correlations. 

Therefore, it is an essential question if the main physical features of high energy electro- 
nuclear reactions can be incorporated in theoretical calculations in an effective manner, 
which would allow one to perform these calculations for a variety of target nuclei in the few 
GeV regime. One essential feature at GeV energies is relativity, and this is the topic with 
which we concern ourselves in this paper. 

Relativistic effects show up in three places: in the kinematics, in the electromagnetic 
current, and in the nuclear/hadronic dynamics. The first point, relativistic kinematics, is 
trivial. In order to treat the nuclear dynamics in a relativistic way, one can solve the Bethe 
Salpeter equation or one of its variations, e.g. the Gross equation 0. The electromagnetic 
current has traditionally been used in a non-relativistic reduction, which assumes that the 
transferred momentum, transferred energy and initial nucleon momentum are small com- 
pared to the nucleon mass. These assumptions are not valid for modern day experiments 
with GeV electron beams and energy and momentum transfers of 1 GeV or more, and there 
exist several improved schemes. In a series of recent publications a new approach to 

the problem of the non-relativistic reduction of the current was introduced, and in , it was 
conjectured that the bulk of the relativistic effects at energies of a few GeV stem from the 
current, not from the nuclear dynamics. This conjecture is based on the large enhancement 
of the transverse-longitudinal response observed in this scheme, which was also seen in the 
fully relativistic calculation of [Q. If this conjecture is indeed true, relativity could be taken 
into account by applying the operator introduced in together with non-relativistic 

wave functions to the calculation of electro-nuclear reactions. This would obviously be very 
useful. 

Although strong evidence was presented in [|l|, this conjecture can only be confirmed 
by a direct comparison to a fully relativistic, microscopic calculation. This comparison is 
carried out in the following sections for the deuteron target. As we wish to concentrate for 
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the moment on the problem of relativity only, we perform this comparison for Plane Wave 
Impulse Approximation (PWIA). A realistic calculation must treat final state interactions, 
too, but the conjecture can be tested successfully at the plane wave level. Moreover, a fully 
microscopic, consistent calculation of final state interaction at GeV energies is extremely 
difficult and not available at the moment. The high energy final state interactions are 
usually treated within the framework of Glauber theory 0. 



II. BRIEF REVIEW OF FORMALISM AND NOTATION 

We start by introducing some notation and giving a brief summary of the basic formalism 
of (e, e'p) reactions. More details can be found in P,|lO . 
The differential cross section in the lab frame is 

( do^ \' _ miMmfPN 

(^Mott free 



\de'dVtf,dVti^ ) Sir^m 

+h{vT'R';l+VTL'RW)], (1) 

where m,, m^r and rrif are the masses of the target nucleus, the ejectile nucleon and the 
residual system, and fl^ are the momentum and solid angle of the ejectile, e' is the energy 
of the detected electron and Qe is its solid angle. The helicity of the electron is denoted by 
h. The Mott cross section is 



2e sm\ej2)^ 



and the recoil factor is given by 



J rec — M- + . [6) 

rrii 

The coefficients vk are the leptonic coefficients, and the Rk are the response functions which 
are defined by 

Rfi ^ \piq)fi\' 

Rl^\Mq}M' + \J4q}M' 
Rf^2^[Jl{q)f,J^{q)f,] 

Ry,^\Mq)M'-\J4Q)M' 

Rf ^ -2^[p\q)f. {.h{q)f. + J-{q)f.)] , (4) 

where the J± are the spherical components of the current. For our calculations, we have 
chosen the following kinematic conditions: the z-axis is parallel to g, the missing momentum 
is defined as pm = Q — PN, so that in Plane Wave Impulse Approximation (PWIA), the 
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missing momentum is equal to the negative initial momentum of the struck nucleon in the 
nucleus, Pm = —p- We denote the angle between and q by 9, and the term "parallel 
kinematics" indicates 6 = 0°, "perpendicular kinematics" indicates 6 = 90°, and "anti- 
parallel kinematics" indicates 6 = 180°. Note that both this definition of the missing 
momentum and the definition with the other sign are used in the literature. In this paper, 
we assume that the experimental conditions are such that either the kinetic energy of the 
outgoing nucleon and the angles of the missing momentum, 6, and the azimuthal angle 0, 
are fixed, or that the transferred energy u, the transferred momentum q, and the azimuthal 
angle 0, are fixed. In the former case, the transferred energy and momentum change for 
changing missing momentum, in the latter situation, the kinetic energy and polar angle of 
the outgoing proton change for changing missing momentum. 
The electromagnetic current operator 



J^(PA;P'A') =u{P'A' 



2m 



N 



u{PA) 



(5) 



where P, P' indicate the four-momenta of the nucleon, can be rewritten in a form that is 
more suitable for application to nuclear problems: 



J^{PA-P'K')^xl J'\P;P') Xa 



(6) 



with 



J' 



p 




M ^1 (^x a) + ^2 (q- ^) (^"'x p) + ^3 [(g X P) ■ a] 



p- 




(7) 
(8) 



Here, fo, ^j' are all functions of u, q^p^; their explicit forms are given in For the reasons 
explained in [||], we refer to the operator associated with as zeroth-order charge operator, 
we call the term containing the first-order spin-orbit operator, the term containing i^i 
first-order convection current, the term containing ^[ zeroth-order magnetization current, 
the term containing first-order convective spin-orbit term, and the term containing ^3 
second-order convective spin-orbit term. 

Note that there are more terms in the full current than with the commonly used strict 
non-relativistic reduction, which assumes that the transferred momentum q is smaller than 
the nucleon mass, and that both the initial nucleon momentum and the transferred energy are 
smaller than q and therefore much smaller than the nucleon mass. Under these assumptions, 
the current operator simplifies to the form 



^nonrel E 



J, 



nonrel 



2mN 



Gm (g X a) + — Ge \P 
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which contains only the zeroth-order charge operator, the zeroth-order magnetization current 
and the first-order convection current. 

In non-relativistic Plane Wave Impulse Approximation (PWIA), we have 



= (^eN free ^(p) > (9) 



where n{p) is the momentum distribution, and the eN cross section is given by 

Erysingle nucleon /■tr\\ 
VkRk (10) 

k 

and the single nucleon responses are related to the nuclear responses by 

j^nudeus ^ (27r)^ j^^ngle nucleon ^^^^ 

SO that one has in total: 

d (7 niN f^f Pn r-i / sr-^ risingle nucleon 

de'dn^dnr, = "(^^ ? "'^^^ 

The momentum distribution is simply the Fourier transform of the wave function, and in 
the non-relativistic case takes the form: 

n{p) = ^{u{pf + w{pf) (13) 

where u{p),w(p) are the S-wave and D-wave functions in momentum space, and the nor- 
malization condition is 

d^pn{p) = 1 (14) 

Note that the different parts of the wave function do not interfere as long as the target 
deuteron is unpolarized. The above formulas, and especially the factorization into eN cross 
section and momentum distribution, is valid only in non-relativistic PWIA, not in PWBA. 
The Born approximation for (e, e'p) scattering contains in addition to the graph where the 
photon couples to to the proton the graph where the photon couples to the neutron. The 
difference between PWIA and PWBA for high energy transfers is small, as the second contri- 
bution involves the deuteron wave function taken at the momentum of the outgoing proton, 
which is high, and therefore the value of the deuteron wave function is very small. However, 
the PWBA is sufficient to break factorization and to permit the interference of the differ- 
ent wave function components. In the following, we present explicit analytic formulas for 
PWIA; these formulas are considerably simpler than the corresponding PWBA expressions. 
For consistency, the numerical results are also presented in PWIA. Their deviation from the 
PWBA results at the kinematics used here is very small, on the order of a few percent. None 
of our results and conclusions are affected by this. 

For the full current operator, the single-nucleon responses take the following form: 
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= ^{Gl + 6'W2) 
= 2Wi + 6^W2 



Rtt = foi^'Sn + - m - ^'Sri) cos(2yp) 

= -6^W2 cos(2^) 

= 2y2cos(v9)-^v/rT7T525iy2 (15) 



with the abbreviations Wi = rGj^ and W2 = i^{G'e + ''"^m)- 

For the strict non-relativistic reduction, the single nucleon responses are given by: 

T3sn 

— 

= 2k^ Gl^ + 6^Gl 
E^^ = -S^Gl cos(2(/?) 

R'^^ = 5G|2\^cos((^) (16) 
The dimensionless variables are defined as follows: 





k1 


K = 




2m N 


6 = 


Pi. 




niN 


T = 




A = 


ijj 






2m AT 



(17) 

The above formulas for the non-relativistic PWIA do not hold for a relativistic treatment 
of the situation. Besides the presence of singlet and triplet P waves, which appear only in a 
full relativistic treatment of the hadronic state, the factorization breaks down. This fact was 



observed earlier, e.g. in ||Tl|--|13|, and can be attributed to the presence of negative energy 
states in the relativistic treatment. To illustrate this point, we write out the expression for 
the hadronic tensor W^^ in the D{e,e'p)n reaction in the relativistic plane wave impulse 
approximation (RPWIA). The relevant kinematic quantities are identified in Fig. |I|; in the 
following, p indicates a four-vector with on-shell energy momentum relation, i.e. p° = Ep := 
y/p'^ + m^: 

W^'' = ^TT{J^{~q)SF{pi + q)r{q) 

XdX2 

Hpu K)KkxM) + ^(-^^1' ^'i)*a;a;a,(Pi)]) (18) 
where the wave function components are given by 
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FIG. 1. This figure shows the relevant kinematic quantities for the electrodisintegration of the 
deuteron in PWIA. 



{'\/2u{p) - w(p))(5ai,--a2] 

V OTT 

^A.A.A.(P) = -^^1:A.-A.(^, ^, -^)2\^[^Piv,{j>)KM + V^^t(p)5A„-Aj , (19) 



where denotes the J = 1 Wigner D function, and the full wave function is 

^A„A,(P) = EN(Pl' Al)vI/tA.A,(Pl) + ^(-Pl, Al)vl/;:^,,,,(pi)] . (20) 
Ai 

Note that whenever the P wave contributions Vg and Vt present in appear in the hadronic 
tensor eq. (0), they are multiplied by a negative energy spinor v. 
After some algebra, the hadronic tensor can be rewritten as 

^ ^Tr(J'^(-g)A+(pi + q)J\q)N) , (21) 

where A"*" is the positive energy projection operator and is the sum of a time-like vector, 
a space-like vector, and a scalar component: 

_ 

with the three distributions given by: 

utv (p) = (p) + (p) + vl (p) + vl (p) (23) 
nsv{p) = u^{p) + w'^{p) - vl{p) - v1{p) 

[{u{p) + ^{2)w{p))vs{p) - {\l{2)u{p) - w{p))vt{p)] (24) 



^P 

ns{p) = u^{p) + w'^{p) - vl{p) - vl{p) 



+ 



2p 
7^ 



m 



[{u{p) + ^{2)w{p))vM - (/(2)^i(p) - w{p))vt{p)] . (25) 



10" 




10"^ - 
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FIG. 2. We show the different momentum distributions as defined in eq. (|25|). The follow- 
ing distributions are presented: the time-like vector distribution (solid line), the space-like vector 
distribution (dash-dotted line), the scalar distribution (dashed line), and the time-like vector dis- 
tribution with the P-waves put to zero (dash-double-dotted line). 



The four vectors used here are the deuteron four- momentum P and k = {0,p), where p is 
the three momentum component of pi. From these formulas it is obvious that apart from 
the time-like vector distribution, there is indeed interference between the P waves and the S 
and D waves in RPWIA. As already noted above, this interference stems from the presence 
of the negative energy terms in eq. (|TB|), and is therefore a genuinely relativistic effect. The 
normalization condition, 

2^A„A^ = / (ifl^ E ^A.A',(p-)r^A.A.(p-) (26) 

leads to the condition: 

1 = - / 7^-Tr(rAr). (27) 
3 7 (27r)3Ep ^ ^ ^ 

Note that the normalization condition in eq. (^) is a relativistic normalization condition 
and differs from the non- relativistic normalization condition in eq. (p^. These differences do 



not matter as long as consistent phase space factors are taken into account in the calculation 
of the cross section. It is clear that only the time-like vector momentum distribution term 
will give a non-vanishing contribution to the trace in the integrand, so that we obtain the 
familiar normalization condition: 

From the hadronic tensor, one can compute the response functions by picking the appro- 
priate values of the indices /i, z/ and evaluating the traces. The explicit expressions for the 
responses are given in appendix In this process, not only the time-like vector momentum 
distribution, but also the scalar and the space-like momentum distributions enter in the eval- 
uation of the response functions. In order to give an idea about the amount of interference 
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FIG. 3. We plot the ratio R = "^^^'^ . The curves indicate the same distributions as in Fig. |2|. 

that is possible in the response functions, we have plotted the three different momentum dis- 
tributions in Figs. ^ and|^. In Fig. ^, we show the time-like vector distribution (solid line), 
the space- like vector distribution (dash-dotted line), and the scalar distribution (dashed 
line). In addition, we have plotted the time- like vector distribution with the P-waves put 
to zero (dash-double-dotted line). This would correspond to the non-relativistic momentum 
distribution. In Fig. |3|, we show the ratio R, i? = , of the deviations of the different 

distributions from the time-like distribution to the time-like vector distribution. One can see 
clearly that both the space-like vector and scalar distributions deviate significantly from the 
time-like vector distribution, in the case of the space-like vector distribution, this happens 
already for rather small momenta. For increasing momenta, the interference effects grow 
considerably, until they reach 30 % to 50%. In contrast, the effect of the P waves in the 
time-like vector distribution is small, and increases only slightly for very high momenta, 
p > 800 MeV/c. This is as expected, as the P waves do not interfere with S or D waves 
in this case, and their normalization is too small to make them significant by themselves, 
compared to the S wave or D wave contribution. 



III. COMPARISON OF DIFFERENT APPROACHES TO RELATIVITY 

In this section, we compare the results of a plane wave calculation of the D{e,e'p)n 
reaction in the framework of the Gross equation (Model R) with a calculation using 
a non-relativistic nucleon-nucleon potential, namely the Argonne V18 potential [|14|, and 



the full relativistic current operator (Model Al). For an overview of all employed models 
and naming conventions, see appendix ^ Both models give a very good fit to the NN 
scattering data. At the high energy and momentum transfer we consider here, the Plane 
Wave Impulse Approximation and the Plane Wave Born Approximation become almost 
identical, they differ by 2% or 3% at most. In order to be consistent with the analytic 
formulas presented in the previous section and in the appendix, the calculations we show 
are done in PWIA. We perform our comparison for the basic observable, the differential 
cross section. 

In Figs. ^ and |^, we compare the fully relativistic, covariant results (Model R) with the 
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FIG. 4. The differential cross section for the ^H{e, e'p)n reaction is shown for fixed 3-momentum 
transfer q = 1.4 GeV/c and different values of the fixed energy transfers lo: 0.61 GeV (a), 0.75 GeV 
(b), and 0.94 GeV (c). The kinematic conditions in the middle panels correspond to quasi- free 
conditions. The solid line shows the result of the fully relativistic calculation (Model R), the dashed 
line shows the result of the calculation with the Argonne V18 wave function and the full relativistic 
current operator (Model Al), and the dash-dotted line shows the result of the calculation with the 
Argonne V18 wave function and the strict non-relativistic reduction of the current operator (Model 
NR). 
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results for the Argonne V18 wave function in conjunction with the full relativistic current 
operator (Model Al). We have also added a curve showing the results obtained with the 
Argonne V18 wave function and the traditional, strictly non-relativistic current operator 
(Model NR), in order to give an idea of the actual size of relativistic effects in the current. 
In the first figure, we have fixed the transferred momentum, q, to q = 1.4 GeV/c, and the 
transferred energy, u, to a value corresponding to the quasi-elastic peak, u = 0.75 GeV 
(middle panel), to a value below the quasi-elastic peak, u = 0.61 GeV (top panel), and to a 
value above the quasi-elastic peak, u = 0.94 GeV (bottom panel). This choice corresponds 
to values of the y-scaling variable ofy = GeV/c, y = —0.2 GeV/c, and y = 0.2 GeV/c. The 



y variable is defined as the negative minimal missing momentum, for details, see e.g. [jT5 
The solid line shows the fully relativistic calculation, the dashed line shows the calculation 
with the relativistic current operator and the AV18 wave function, and the dash-dotted line 
shows the the calculation with the traditional non-relativistic reduction of the current and 
the AVIS wave function. We start our discussion with the lowest energy transfer, u = 0.61 
GeV. One sees that the fully relativistic calculation and the AVI 8 with full relativistic 
current are relatively close to each other, with the fully relativistic result a bit lower for 
missing momenta around 400 MeV/c. The result obtained with the non-relativistic current 
operator is larger than the relativistic curves. For the quasi- free kinematics (middle panel), 
the fully relativistic calculation and the AVIS plus full relativistic current agree quite nicely, 
with the fully relativistic curve being slightly lower for momenta around 350 MeV/c and 
then being slightly higher for momenta higher than 500 MeV/c. The non-relativistic result 
is significantly higher for all missing momenta. This trend continues for the kinematics 
above the quasi-elastic peak (lower panel), where the difference between non-relativistic and 
relativistic curves increases again. This is to be expected, as one of the basic assumptions of 
the strictly non-relativistic reduction is that the transferred energy is much smaller than the 
transferred momentum. This condition is fulfilled to a certain extent below the quasi-elastic 
peak, but not on or above it. The fully relativistic result and the AVIS plus full relativistic 
current agree for missing momenta up to 400 MeV/c, then the AVIS result drops off a bit. 
This difference is the largest between these two calculations which we have encountered 
so far, and it is logical that it should arise where it does: certainly, relativistic effects are 
expected to be strongest for situations where both energy and momentum transfer are large, 
and for high missing momenta, which correspond in the absence of final state interactions 
to the initial momentum of the nucleon in the nucleus. It is logical to expect any relativistic 
effect in the nuclear dynamics to show up at high missing momenta. Below, we will discuss 
the relativistic effects surpassing kinematics and current. 

In order to get a different view, we also show the same results for a different kinematic sit- 
uation: now, we fix the kinetic energy of the outgoing proton to 1 GeV, and consider parallel 
and perpendicular kinematics. The results are shown in Fig. For parallel kinematics (top 
panel), the fully relativistic result and the AVIS plus full relativistic current result coincide 
for missing momenta below 300 MeV/c, for higher missing momenta, the fully relativistic 
curve is a bit lower. The non-relativistic result differs from the other results already at the 
lowest missing momenta, but the difference is comparatively small in these kinematics, the 
momentum transfer increases very quickly, so that the relation u << q is roughly applicable, 
as in the case of Fig. ^, top panel. In perpendicular kinematics, the fully relativistic result 
and the AVIS plus exact current result agree very nicely up to Pm = 450 MeV/c, for higher 
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FIG. 5. The differential cross section for the '^H{e, e'p)n reaction is shown for fixed kinetic 
energy, T^in = 1 GeV, of the outgoing proton in parahel and perpendicular kinematics. The curves 
have the same meaning as in Fig. 

missing momenta, the two calculations start to diverge, leading to a difference similar to the 
one observed for the highest energy transfer in the fixed q kinematics. 

To summarize the results of our comparison so far, we have seen that there is very good 
overall agreement between the fully relativistic result (Model R) and the AV18 plus full rel- 
ativistic current result (Model Al), with some slight deviations at higher missing momenta, 
Pm > 450 MeV/c, especially for energy transfers u comparable to q. The calculation with 
the AV18 wave function and the non-relativistic reduction of the current (Model NR) differs 
from the relativistic results considerably, starting from the lowest missing momenta. In the 
light of these results, it is fair to say that the bulk of the relativistic effects in the few GeV 
region does stem from the current, as conjectured in 

In order to ascertain the influence of relativity in contrast to the influence of the modeling 
of the NN interaction properly, we carry out an additional comparison. We take a parame- 
terization of the S wave and D wave part of the Gross wave function, and use this together 
with the full relativistic current (Model Gl). Whether the model for the A^A^ interaction 
is relativistic or not, this choice alone does not specify the model. The data constrain the 
deuteron wave function somewhat, but it is well known that for higher momenta, where the 



D wave is dominant, different models make different predictions. E.g. the Bonn [16| and 



the Paris |T^ wave functions are both non-relativistic, but differ greatly in their D wave 
content. In our comparison, we are interested in the different treatment of relativity, and 
want to eliminate differences due to the difference in the A^A^ interaction modeling. We 
therefore show a plot of the momentum distribution as given by eq. (|13D of several non- 
relativistic models, and also plot the momentum distribution obtained with the S wave and 
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p (MeV/c) 



FIG. 6. The momentum distribution n(p) for the deuteron is shown for several wave func- 
tion models. The solid line shows the Bonn result, the dashed line shows the Paris result, the 
dash-dotted line shows the CD Bonn result, and the long dashed line shows the Argonne V18 
result. The short-dashed line shows the momentum distribution calculated with the S wave and D 
wave of the Gross wave function. 

D wave of the Gross wave function. We omit the P waves here in order to compare the same 
quantity, the different normahzation induced by this has been taken into account. As seen 
in Fig. all wave function models agree for momenta up to 300 MeV/c, and then start 
to diverge. The two Bonn potential curves, Bonn and CD Bonn |jl8[, are lower than 



Paris, Argonne V18 and the Gross wave function. The latter three are fairly close, with the 
Paris and Argonne V18 slightly higher than the Gross wave function. This means that our 
comparison between the fully relativistic calculation (Model R) and the Argonne V18 with 
the full relativistic current (Model Al) is meaningful, as the differences in these calculations 
should arise mainly from the different treatment of relativity, not from the differences in the 
wave function. Note that a comparison of the fully relativistic calculation with one of the 
Bonn potential wave functions plus full relativistic current would not have achieved that. 

In order to be as precise as possible in our investigation of relativistic effects, we compare 
the fully relativistic calculation (Model R), the Argonne VIS plus full relativistic current 
(Model Al), and the Gross S wave and D wave plus full relativistic current (Model Gl). The 
curves are shown in Figs. |^ and ^. The kinematic conditions are the same as in Figs. ^ and 
^. First, we discuss the results for fixed energy and momentum transfer. Due to the fact 
that the momentum distribution calculated with the Gross wave function is slightly lower 
than the Argonne V18 momentum distribution, the Gross S+D plus full relativistic current 
results are slightly lower than the AV18 plus full relativistic current results. This improves 
the agreement with the full relativistic calculation for the lowest energy transfer, uj = 0.61 
GeV for momenta less than 400 MeV/c. For higher missing momenta, the deviation from 
the full relativistic calculations increases a bit, although not significantly, leaving a little 
more room for genuine relativistic nuclear dynamics effects. In Fig. |], we show the results 
for fixed kinetic energy of the outgoing proton. For parallel kinematics, the use of the Gross 
S and D waves plus full relativistic current results improves the agreement with the fully 
relativistic calculation, the two curves differ only slightly. For perpendicular kinematics, all 
three curves agree nicely for missing momenta up to 350 MeV/c, then the Gross S and D 
wave plus full relativistic current calculation differs from the full calculation a bit more than 
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FIG. 7. The differential cross section for the ^H[e, e'p)n reaction is shown for fixed 3-momentum 
transfer q = 1.4 GeV/c and different values of the fixed energy transfers u: 0.61 GeV (a), 0.75 GeV 
(b), and 0.94 GeV (c). The kinematic conditions in the middle panels correspond to quasi- free 
conditions. The solid line shows the result of the fully relativistic calculation (Model R), the dashed 
line shows the result of the calculation with the Argonne V18 wave function and the full relativistic 
current operator (Model Al), and the dash-dotted line shows the result of the calculation with the 
Gross S wave and D wave plus the full relativistic current operator (Model Gl). 
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FIG. 8. The differential cross section for the '^H{e, e'p)n reaction is shown for fixed kinetic 
energy, T^in = 1 GeV, of the outgoing proton in parahel and perpendicular kinematics. The curves 
have the same meaning as in Fig. ^. 

the corresponding AV18 curve. 

So, overall, we have very good agreement between the fully relativistic calculation and the 
AV18 or Gross S and D wave plus full relativistic current. In contrast to this, a calculation 
with the strict non-relativistic reduction of the electromagnetic current disagrees markedly 
with the relativistic results. In the remainder of this paper, we discuss the possible sources 
for the remaining small discrepancies at higher missing momenta. 



A. P-wave contributions 

One genuinely relativistic effect that cannot be taken care of by the electromagnetic 
current is the emergence of P waves in the deuteron wave function. They arise naturally from 
the lower spinor components, and there is a P-wave singlet and a P-wave triplet contribution, 
see e.g. |]19|. Although their normalization is extremely small compared to the S-wave and 



the D-wave, they can play an important role, e.g. in the calculation of the elastic deuteron 
form-factors, where they shift the location of the minimum of B^Q"^) by a considerable 
amount . How important are they for (e, e'p) reactions? We investigate this question 
by switching off the P-wave contribution and comparing the full result (Model R) with the 
result without P-waves (Model RSD) in Figs. ^ and Overall, the differences between 
the full calculation and the calculation without P waves are very small. In the fixed q and 
uj kinematics, the two curves coincide for the lowest energy transfer of 0.61 GeV, and barely 
differ for uj = 0.75 GeV. Only for the highest energy transfer of 0.94 GeV, the P waves have 
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an effect, their inclusion slightly increases the result. We would like to point out that this 
is due to interference between the P waves and the S and D waves. The contribution of 
the P waves alone is approximately two orders of magnitude smaller than the contribution 
of S and D waves, so the interference effects must be quite large in order to be visible. In 
Fig. in parallel kinematics (top panel), the two curves barely differ. It is interesting, 
though, to see that the interference in this case is destructive, leading to a slightly enhanced 
result when the P waves are omitted. In perpendicular kinematics, the omission of the P 
waves lowers the result slightly for missing momenta higher than 400 MeV/c. In each case, 
when comparing with Figs. ^ and ^ and with Figs. and |^, one sees that the P wave 
contribution explains a part of the difference between the fully relativistic result and the 
other calculations, but not all of it. Another important aspect for calculations at higher 
missing momenta are off-shell prescriptions, they will be discussed in the next section. 

One additional remark is in order: although the P-waves are only of marginal importance 
for the calculation of the differential cross section and the large responses, they may play a 
more important role in the calculation of single or double polarization observables, by inter- 
fering with one of the larger contributions. Similar situations arise e.g. when considering the 
second-order convective spin-orbit contribution to the transverse-transverse response func- 
tion, which is considerable due to the interference with the magnetization current, the largest 
current component that cannot contribute to the transverse-transverse response function 
otherwise ^] or when considering the contribution of the spin-orbit final state interaction 
to the fifth response function, which allows for a large interference of the two dominant com- 
ponents of the current that is not present otherwise 0]. In general, one needs to be careful 
in situations where a change of the spin structure will allow for a new, large contribution. 
In these situations, it is advisable to take a close look at the relevant Clebsch-Gordan coeffi- 
cients and to determine if additional important contributions could arise from the presence 
of P-waves. 



B. OfF-shell effects 

Whenever one applies the one-body electromagnetic current operator not to a free nu- 
cleon, but to a nucleon bound in the nucleus, one needs to introduce an off-shell prescription. 
Nucleons in the nucleus are bound and therefore off-shell; they do not fulfill the same energy- 
momentum relations as free nucleons. Currently, there exists no microscopic description of 
this off-shell behavior that can be applied for a wide range of kinematic conditions — there 



are only ad hoc prescriptions, which lead to differing results for certain kinematics pl|-p3 
The variations tend to increase with increasing momentum of the initial nucleon - the higher 
its momentum, the further it is off-shell. The fully relativistic calculation (Model R, Model 
RSD) uses a more sophisticated off-shell prescription, which takes into account the actual 
kinematics. In the calculation with the non-relativistic wave functions and the exact form of 
the current operator, we have chosen the popular ansatz of employing the "on-shell form" of 
the current. However, these are specific choices among an infinite number of prescriptions. 
We will quantify the theoretical uncertainty due to these choices in the following. It is easy 
to see why there are, in principle, infinitely many possibilities to choose an off-shell prescrip- 
tion. Consider the case of non-relativistic PWIA, and of factorization holding there, i.e. we 
are investigating the off-shell effects in the framework of non-relativistic wave function and 
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FIG. 9. The differential cross section for the ^H{e, e'p)n reaction is shown for fixed 3-momentum 
transfer q = 1.4 GeV/c and different values of the fixed energy transfers u: 0.61 GeV (a), 0.75 GeV 
(b), and 0.94 GcV (c). The kinematic conditions in the middle panels correspond to quasi- free 
conditions. The solid line shows the result of the fully relativistic calculation (Model R), the dashed 
line shows the result of the fully relativistic calculation without the P waves (Model RSD). 
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FIG. 10. The differential cross section for the ^H{e, e'p)n reaction is shown for fixed kinetic 
energy, Tfcj„ = 1 GeV, of the outgoing proton in parallel and perpendicular kinematics. The curves 
have the same meaning as in Fig. ^. 



full current operator now. Due to the factorization, the cross section can be described by 
the product of momentum distribution and off-shell electron-proton cross section, see eq. 
(y). For a free nucleon, one can write down the ep cross section, and then rewrite it using 
simple algebra, that in general will employ the on-shell energy momentum relation for the 
free nucleon. On-shell, all these ways of writing the ep cross section are equivalent - off-shell, 
they differ, and of course there is no way to identify a single expression as the correct one. 
E.g., one obtains different versions of the off-shell ep cross section if one starts out from 
the one-body current in the form given in eq. (^, or if one uses the Gordon decomposition 
on the current. This is precisely how deForest [|21[ obtained his off-shell prescription ccl 
(current like in eq. (|^)) and cc2 (Gordon decomposition applied to the current). At first 
sight, the situation might look rather grim, however, things are not as bad as one might as- 
sume. In the following, we compare the employed off-shell prescriptions. In eq. (^, we gave 
the expressions for the single nucleon responses in terms of the ^s, as they have been used 
in the calculations shown above (Models Al and Gl), and also in terms of Wi, W2 (Model 
A2). This latter form was obtained by doing some algebra involving the on-shell kinematic 
relations, and therefore constitutes a different off-shell prescription. In Fig. 11, top panel, 
we show the ratio of the ep cross section calculated with the ^s, called off-shell prescription 
1, and with Wi,W2, referred to as off-shell prescription 2. The curves show the ratios for 
the kinematic settings discussed in this paper. One sees that, according to expectations, the 
off-shell prescriptions agree for momenta below 200 MeV/c, and then start to diverge. At 
the highest missing momenta considered here, the deviations range from 10% to 25%. At 
the intermediate momenta, the typical deviation from 1 is about 5 %. The largest deviations 
occur for T^j^ = ^GeV in perpendicular kinematics (long dashed line), and for the q = 1.4 
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GeV/c, uj = 0.61 GeV case (dash-dotted line). In the middle panel of Fig. |TT], we show the 
ratio of off-shell prescription 1 and the off-shell prescription used in the covariant calculation, 
denoted as off-shell 3. The explicit formulas for off-shell 3 are given in appendix y. Again, 
the two off-shell prescriptions agree nicely for low missing momenta, and start to diverge for 
Pm > 200 MeV/c. The deviations are large, almost 50 % at the highest missing momenta for 
Tkin = ^GeV in perpendicular kinematics, and also considerable at high missing momenta 
for q = 1.4 GeV/c, u = 0.61 GeV. In all other considered kinematics, the deviations do 
not exceed 5% even at the highest missing momenta. Note that the deviations from 1 for 
the ratio of off-shell 2 to off-shell 3 are smaller than any of the ratios shown here, we have 
chosen to display the "extreme" cases here. 

Although the off-shell prescriptions 1,2, and 3 are all relativistic, there remains an am- 
biguity in the choice of the off-shell kinematics. The scope of this paper was to investigate 
the treatment of relativistic effects, not the specific off-shell choices. However, our results 
give a good idea about the size of the uncertainties, and about the kinematic regions which 
are sensitive to them; whereas in parallel kinematics, the results are quite insensitive to the 
employed off-shell prescription, the results in perpendicular kinematics are more sensitive. 
Below missing momenta of 200 MeV/c, off-shell ambiguities practically do not exist, and for 
Pm < 400 MeV/c, they will in general not be larger than 10 %. 

Sometimes, the full relativistic current is rejected as introducing off-shell contributions, 
and the claim is made that the non-relativistic reduction is preferable as it would avoid such 
ambiguities. This argument is not correct. Firstly, even the non-relativistic version of the 
current contains a dependence on the nucleon momentum, leading to terms proportional to 
the perpendicular momentum squared in the single nucleon responses, see eq. (|16|). On-shell, 
this momentum can be replaced with the energy by using the on-shell kinematic relations, 
and therefore one does have off-shell ambiguities even in this case. One might argue that 
the non-relativistic versions of the charge operator and the magnetization current do not 
contain any dependence on the initial nucleon momentum, and should therefore be free of 
off-shell ambiguities. Here, one needs to ask, however, if this feature was not bought at a 
too high price: the ep cross sections calculated in the strict non-relativistic reduction differ 
from the relativistic version by 30 % to 80 %, and this discrepancy starts at p^ = MeV/c. 
We demonstrate this point by presenting the ratio of the relativistic to the non-relativistic 
e—p cross section for one specific off-shell prescription, namely the on-shell form with the ^s, 
referred to as off-shell 1, see Fig. [Tl|, lower panel. Even though the off-shell effects in the full 
relativistic version of the current do introduce a theoretical uncertainty at higher missing 
momenta, the relativistic treatment is definitely superior to the non-relativistic version, 
where one knows that the results are off for all missing momentum values, and they are off 
up to a factor of 5. 



IV. SUMMARY 

The principal goal of this paper has been to investigate the origin of the relativistic 
effects in the electrodisintegration of the deuteron. Apart from the obviously relativistic 
kinematics, relativistic effects occur in the electromagnetic current operator and the nuclear 
dynamics. We have shown that calculations with a non-relativistic wave function and the 
fully relativistic current operator reproduce the fully relativistic calculation very well up to 
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FIG. 11. The upper panel shows the ratio of the off-shell electron proton cross sections calcu- 
lated with prescription 1 and prescription 2. The different lines correspond to different kinematic 
settings: the solid line shows the ratio for T^in = 1 GeV, parallel kinematics, the long-dashed line 
for Tkin = 1 GeV, perpendicular kinematics, the dash-dotted line for q = 1.4 GeV/c, to = 0.61 
GeV, the dotted line for q = 1.4 GeV/c, co = 0.75 GeV, and the short-dashed line for q = 1.4 
GeV/c, uj = 0.94 GeV. In the middle and the lower panel, the different lines correspond to the 
same kinematic conditions as in the upper panel of this figure. The middle panel shows the ratio 
of the off-shell electron proton cross sections calculated with prescription 1 and prescription 3. 
The lower panel shows the ratio of the off-shell electron proton cross sections calculated with the 
full relativistic current and the strict non-relativistic reduction of the current for a fixed off-shell 
prescription, namely off-shell 1. 
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400 MeV/c of missing momentum in all kinematics. At the same time, a calculation with 
a non-relativistic wave function and a strictly non-relativistic current operator drastically 
fails to reproduce the relativistic results, in all kinematics and even for the lowest missing 
momenta. 

We have shown that the relativistic effects in the nuclear dynamics, namely, the presence 
of P-waves, are rather small in the cross section. They do appear only in places where one 
would expect relativistic effects to show up, i.e. for high missing momenta, pm > 400 MeV/c 
and for energy transfers which are large and comparable to the momentum transfer. The 
remaining discrepancy that we found in a few kinematics for high missing momenta was 
shown to stem from the different off-shell prescriptions used in the two calculations. 

Since shows that the effects of bound state dynamical relativity are small also in 
heavier nuclei, the technique of can reliably be apphed to heavier nuclei, too. 

These findings verify the conjecture of ^ that at energies of a few GeV, the bulk of the 
relativistic effects in electronuclear reactions stems from the current operator, and not from 
the nuclear dynamics. 



ACKNOWLEDGMENTS 

The authors thank J. Adam, T. W. Donnelly, and F. Gross for many stimulating discus- 
sions on the subject of relativity. The authors thank R. Schiavilla for providing a parameter- 
ization of the Argonne V18 and CD Bonn deuteron wave functions. This work was in part 
supported by funds provided by the U.S. Department of Energy (D.O.E.) under cooperative 
research agreement #DE-AC05-84ER40150. 

APPENDIX A: OVERVIEW OVER PRESENTED CALCULATIONS: NOTATION 



For the convenience of the reader, we list all calculations and assign a short name for 
them: 



wave function 


responses/off-shell prescription 


description 


name 


Gross wave function; 
S, D, P waves 


as in appendix 
off-shell 3 


fully relativistic, 
covariant calculation 


Model R 


Gross wave function; 
S, D waves 


as in appendix 
off-shell 3 


fully relativistic, 
covariant calculation 


Model RSD 


Argonne VI 8; 
S, D waves 


as in eq. (15), 
with the ^, off-shell 1 


AV18 plus full 
relativistic current 


Model Al 


Argonne V18; 
S, D waves 


as in eq. (|15|), 
with Wi,W2, off-shell 2 


AV18 plus full 
relativistic current 


Model A2 


Gross wave function; 
S, D waves 


as in eq. (|15|), 
with the ^, off-shell 1 


Gross S and D wave function 
and full relativistic current 


Model Gl 


Argonne V18; 
S, D waves 


as in eq. (|16|), 
with the ^, off-shell 1 


AV18 plus strictly 
non-relativistic current 


Model NR 
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APPENDIX B: FULLY RELATIVISTIC NUCLEAR RESPONSES 

Here, the explicit analytic expressions for the nuclear response functions using the off- 
shell prescription of the fully relativistic, covariant calculation are given: 

= ^ ^ {(CrriTvip) + Crnsvip) + C^ns{p)) 
+Cl''nsv{p) [(M,-£;,)|-|p,]2 

+ {CrnTv{p) + C'/nsv{p)) 2M, | [(M, -E,)^- (Bl) 

Rt^^ {-2 {CruTvip) + Cf'^nsvip) + Cfusip)) 

+Cfnsv{p) pI} (B2) 
Rtt = C^'^nsvip) pI cos(299) (B3) 

OTT 

i?TL = ^ I {Crnsvip) 2V2 [(M, - - P± cos((^) 

+ {CrnMp) + Cf "'n^vb)) 2V2 M, cos(<^)} (B4) 

The coefficients C of the time-like vector, space-hke vector, and scalar momentum dis- 
tributions are given by: 

+2M,cuF,{Q') {F,{Q') + F,{Q')) - ^^^F^{Q')} (B5) 
Cj^ = (B6) 

Cj^ = (B8) 
where A — pf — — M| — 2MaEp measures the "off-shellness" of the bound nucleon. 

+|||j'.W')(f'iW')+-F'2(0')) (B9) 

^ 777' 

= A + (BIO) 
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2 m' 



:F,{Q')F,{Q' 



(Bll) 

(B12) 

(B13) 

(B14) 
(B15) 
(B16) 



Note that the responses given here are to be inserted into eq. (|ID for the cross section in 
the lab frame. The responses R used in |T^ in eq. (95) are connected to these responses in 
the following way: 



Rl 
Rx 
Rtt 
Rtl 



Rl 



27r^ W 
m? MtQ^ 

——R 

27?^ W ~ 

— TT^TT cos(2(/?) 
27r2 



Rtl cos(v9) 



m2 MtQ 

Here, W indicates the invariant mass of the final state, and Mt is the target mass 



(B17) 



APPENDIX C: SINGLE NUCLEON OFF-SHELL RESPONSES 

In the case of vanishing P-wave contributions, the different momentum distributions 
coincide: 

nTv{p) = nsv{p) = ns{p) = n^{p) = u^{p) + (CI) 

Note that the normalization condition of the Gross wave function as given in eq. (^) differs 
from the one employed for the Bonn and Paris wave function, eq. (0). The different 
normalizations have of course been taken into account when we compared the two results. 
For vanishing P-waves, the responses factorize, and the resulting off-shell single nucleon 
responses are called "off-shell 3" in the text. For the convenience of the reader, they are 
given here explicitly: 

RT"" = ^ ^ {-2r(Fi(g2) + F,mf + ,r^,{Fim + rF|(g^)) - ^FKQ^) 



23 



^ {Fl{Q^)+TFi{Q^))2MA[{Md-E,)^-'^p,]} (C2) 



+^,i^iiQ') + rFi{Q'))pl} (C3) 

it^TT"' = + ^^2 cos(2(^) (C4) 
Sir 

^ :{F^{Q') + rF|(g2)) 2V2 M, l-p^} cos(vp) (C5) 



m 
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